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Abstract: There is a shock wave in the supercritical airfoil when the inflow velocity is faster than the 


critical Mach number, which leads to the wave drag increase rapidly. In order to reduce the wave drag for 
different flow fields, the shock control bump are studied in this thesis.An optimization platform was 
developed consisting of a genetic algorithm optimiser, a NURBS curve module and a CFD simulation 
module and the different flow states were optimized by using this platform.A calculation for RAE2822 
airfoil has been done, and the computation results demonstrate that the optimized bump configuration can 
reduce the wave drag of the airfoil greatly and improve the performance for off-design conditions. To 
overcome the problem of narrow effective regions of significant drag reduction for one-point designed 
bumps, a design of two-dimensional adaptive bump based on shape memory alloy is proposed. Shock 
wave in different flow condition can be controlled through adjusting the appearance of SMA bump by 


temperate. 
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部 激 波 出 现 ， 且 激 波 强度 随 着 飞行 速度 的 提高 迅速 


osl 增加 。 根 据 波音 公司 的 研究 , DUR EF CH (BWB) 
J&4 Wr, 2025 年 的 全 球 航空 运输 量 将 比 目 前 增 布局 为 代表 的 下 一 代 民 机 相 比 现 有 常规 机 身 + 机 距 

长 2-3 倍 。 根 据 美 国 “ 新 一 代 航 空运 输 系 统 外 形 ， 在 噪声 、 油 耗 等 方面 有 望 同时 满足 新 一 代 先 

CNextGen)” 计 划 ， 预 计 在 2030-2035 年 投入 使 用 进 民 机 的 技术 要 求 。 研 究 表明 ，BWB 构 型 的 激 波 

的 先进 民 机 一 一 N+3 代 ， 要 求 其 排污 、 油耗 、 运 营 阻力 显著 增 大 ， 大 于 总 阻力 的 10% 中 。 因 此 ， 针 对 

成 本 等 方面 均 有 较 大 幅度 的 改善 趾 。 翼 型 的 激 波 控制 减 阻 技术 有 很 强 的 应 用 前 景 。 
飞行 器 在 高 亚 音速 飞行 状态 , BUR. BLUE A8 I9) 鼓 包 减 阻 技术 , 是 一 种 典型 的 用 于 激 波 控制 的 流 
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动 控制 手段 。 已 有 的 研究 成 果 表 明 ， 实 体 鼓 包 能 够 
显著 弱化 激 波 强度 、 减 小 激 波 阻力 ， 同 时 不 附加 黏 
性 阻力 FE-9。 
国内 外 多 个 航空 研究 单位 均 对 鼓 包 激 波 减 阻 特 
性 进行 了 研究 。1992 年 ，Ashill 等 人 第 一 次 提出 了 
在 层 流 嗣 型 上 使 用 鼓 包 减 小 激 波 阻 力 的 思路 门 。 
NASA 兰 利 研 究 中 心 针 对 高 亚 音速 条 件 ， 对 安装 了 
固定 式 鼓 包 的 辟 型 在 非 设计 点 时 的 气动 状态 进行 了 
风 洞 实验 研究 ,结果 表明 通过 选用 合理 的 鼓 包 构 型 ， 
翼 型 阻力 能 够 减 小 1296-1590, 英国 Glasgow 大 学 
K.H..LO 等 人 在 超 音 速 (Mach=1.3) 来 流 条 件 下 ， 
研究 了 三 维 鼓 包 构 型 在 经 型 流动 分 离 控 制 方面 的 应 
。 结 果 表 明 ， 三 维 鼓 包 带 来 的 喷射 压力 有 效 的 控 
制 了 气流 方向 和 顺 咽 展 方向 上 的 流动 分 离 趾 。 剑 桥 
大 学 Jereny P.Eastwood 等 人 对 三 维 辟 包 构 型 设计 进 
行 了 研究 ， 分 析 了 不 同 设计 参数 对 鼓 包 减 阻 效 果 的 
影响 5 。 


鼓 包 减 阻 的 多 项 研究 表明 ,固定 式 鼓 包工 作 带 宽 


= 


宕 ， 难 以 实际 应 用 。 为 此 研究 人 员 提 出 了 自 适 应 鼓 
包 的 设计 概念 ， 希望 根据 不 同 的 飞行 状态 自动 调整 
喜 包 构 型 ， 以 满足 不 同 条 件 下 的 飞行 要 求 。E.Jinks 
等 人 研究 了 单 点 驱动 的 铝 制 自 适应 豆包 结构 ， 结 果 
表明 自 适 应 鼓 包 能 够 对 多 个 任务 状态 进行 优化 ， 可 
保证 器 型 在 更 大 的 包 线 内 具有 较 优 的 气动 特性 1。 
加 拿 大 AV Popov, M Labib 等 人 研究 了 分 布 式 驱 动 
的 自 适 应 鼓 包 变 形 控制 ， 对 鼓 包 变 形 过 程 中 驱动 器 
控制 方法 、 变 形 驱 动力 等 进行 了 分 析 呈 。 

本 文 以 超 临 界 缀 型 激 波 减 阻 为 目标 ， 结 合 
NURBS 曲线 建 模 、CFD 仿真 等 方法 研究 以 
RAE2822 为 代表 的 超 临 界 辟 型 鼓 包 位 置 、 不 同 流动 
条 件 下 的 喜 包 构 型 的 优化 问题 ， 并 在 此 基础 上 提出 
了 一 种 自 适 应 鼓 包 结构 实现 形式 。 


1 鼓 包 构 型 设计 


1.1 NURBS 曲线 

为 了 得 到 鼓 包 在 不 同 气动 状态 下 的 外 形 ， 利 用 
样 条 曲线 对 鼓 包 进 行 参数 建 模 。 常 用 的 曲线 参数 设 
计 方 法 包括 : Bezier 曲线 、B 样 条 、NURBS 
(Non-Uniform Rational B-Splines， 非 均匀 有 理 B 样 
条 ) 等 由， 可 以 根据 控制 点 来 生成 目标 几何 模型 。 
HT NURBS 曲线 具有 良好 的 局 部 支撑 性 ， 即 改变 
线 某 一 控制 点 的 值 只 会 对 与 之 对 应 的 曲线 局 部 产 
生 影响 ， 从 而 实现 了 对 几何 模型 的 局 部 修改 "4。 本 
文选 用 NURBS 曲线 作为 鼓 包 参数 化 方法 ， 利 用 
MATLAB 编程 实现 自 适应 鼓 包 构 型 的 参数 化 建 模 。 
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通常 ， 用 参数 u 表示 的 p 次 NURBS 曲线 点 定 
义 为 
PN, (00) P æy) 


C(u) = += , a<u<b (1) 


其 中 ，{P yy 是 控制 点 (构成 曲线 的 控制 多 


边 形 )，ww, 表 示 权 重 因 子 ， 用 于 定义 控制 点 的 权重 。 


nl 个 权 因子 分 别 于 同等 数量 的 控制 点 (P. Gy) 


相对 应 。 {N,, 00))} 是 由 节点 矢量 UV ， 根 据 递 推 关系 
得 到 的 p 次 B 样 条 基 函 数 。 节 点 矢量 
U = 也 ,22 由 一 组 实数 序列 构成 ， 习 为 节点 


RE, u 为 节点 值 ， 


文选 用 均匀 节点 矢量 进行 描述 。 


u, uj, 120,1, m-l. A 


N, (u) - i i+l 
0 
N, (u)= (2) 
hp 
u— u, Ui p TU 
N, sae) —— —N a (u) 
itp ^ "i 计 P+L — "ik 
U —1d,-- d, usur uus By b] 
pil pH (3) 
Ui Eu 


因此 ， 人 需 给 出 组 控制 点 {PG%,»)) 的 坐标 


寺 ， 指 定 曲线 的 阶 次 、 各 个 控制 点 的 权重 ， 就 可 以 
得 到 一 条 NURBS 曲线 。 如 图 1 所 示 ， 为 利用 4 个 
控制 点 生成 的 3 阶 NURBS 曲线 ， 各 控制 点 的 权重 
均 为 1。 图 中 ， 曲 线 01 控制 点 2 的 Y 坐标 为 0.5, 
hek 02 控制 点 2 的 Y 坐标 值 为 1.0。 可 以 看 出 ， 通 
过 改变 控制 点 的 坐标 值 ， 可 以 对 曲线 形状 进行 局 部 
调整 。 


L5 一 一 一 NURBS-01 Y2-0.5 
一 @@ 一 控制 点 Y12=0.5 


= = = NURBS-02 Y2=1.0 
一 重 一 控制 点 Y22=1.0 


图 1 NURBS 曲线 随 控制 点 的 变化 

Fig. 1 NURBS Curve with the change of control points 
1.2 激 波 控制 鼓 包 构 型 优化 方法 

为 了 得 到 特定 流 场 条 件 下 的 鼓 包 构 型 ， 本 文 使 
遗传 算法 对 鼓 包 NURBS 曲线 参数 进行 优化 。 利 
] CFD 仿真 得 到 翼 型 的 气动 特性 参数 〈 升 力 系数 、 
阻力 系数 、 升 阻 比 等 )。 在 保证 升力 系数 的 条 件 下 ， 
寻求 具有 最 小 阻力 系数 的 鼓 包 构 型 。 
利 包 曲线 生成 模块 ， Gambit 网 格 更 新 模块 ， 
Fluent 翼 型 气动 参数 计算 模 用 isight 多 学 科 仿真 优 
化 平台 ， 通 过 集成 Matlab 鼓 块 ， 使 用 GA IDEE 
化 算法 ,实现 了 对 不 同 流动 状态 下 的 鼓 包 构 型 优化 。 
流程 图 如 图 2 所 示 


Optimizer 
ISIGHT E 
目标 方程 
GA 
NURBSÉ zE 曲线 
设计 变量 
Output File 
Ci, Ca LD ... 
New Airfoil 
更 新 CFD 计 算 网 格 流 场 求解 
Gambit Fluent 


DS 


2 SRCRLIAEEL (LAGU: 
Fig.2 Flow chart of optimization of bump configuration 
1.3 CFD 可 靠 性 验证 
CFD 仿真 技术 广泛 应 用 于 多 种 主动 、 被 动 流动 
控制 技术 的 研究 工作 中 ， 本 文选 用 Fluent 流体 仿真 
软件 对 RAE2822 辟 型 鼓 包 减 阻 的 气动 特性 进行 研 
究 。 为 方便 起 见 ， 选 用 RAE2822 SUV EC TERES 
subire 能 力 进行 验证 。 
计算 状态 : 马赫 速 Ma-0.729, I AOA=2.31° , 


雷诺 数 Re=6.49981e+06。 

K 1 为 Fluent 计算 结果 与 实验 结果 的 对 比 。 可 
以 看 出 ， 采 用 S-A 模型 计算 的 阻力 偏 大 ，SST 模型 
计算 的 升力 偏 小 。 由 于 优化 过 程 需 要 迭代 多 次 ， 综 
合 考虑 计算 时 间 、 计 算 精 度 。 选 用 S-A 模型 进行 后 
续 计 算 ， 网 格 数 59800， 第 一 层 网 格 高 度 6.3X10” 
倍 弦 长 。 

表 1 RAE2822 升 阻 特性 计算 与 试验 结果 对 比 

Table 1 Comparison of calculation and experimental 


results for RAE2822 aerofoil 


Grid number Cd Ci 
Experiment 0.0121165 0.730962 
59800 0.0138742 — 0.700442 
S-A 64400 0.0138196 — 0.700179 
75600 0.0138010 — 0.699636 
59800 0.0132308 0.677409 
k-o SST 64400 0.0131077 — 0.678579 
75600 0.0130952 — 0.678749 


如 下 图 所 示 为 S-A 模型 计算 得 到 的 压力 系数 分 
布 ， 结 果 表 明 在 上 述 条 件 下 计算 得 到 的 RAE2822 
辟 型 表明 压力 系数 分 布 与 NASA 的 风 洞 试验 测试 结 
果 吻 合 较 好 。 从 计算 结果 可 以 看 出 ， 本 文采 用 的 流 
场 计算 方法 是 可 靠 的 。 
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Fig.3 Comparison of pressure coefficient between results of 


"s 


simulation and experiment 
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不 同 的 鼓 包 变形 区 域 、 变 形 形 状 会 带 来 不 同 的 
气动 效果 。 为 了 确定 合理 的 局 部 鼓 包 变形 区 域 ， 在 
辟 型 上 表面 激 波 位 置 附近 ， 选 取 三 个 不 同 的 优化 区 
间 用 于 确定 合理 的 鼓 包 变 形 区 域 ， 分 别 为 Case 1, 
Case 2, Case 3， 如 图 4 所 示 。 变 形 区 间 取 值 范围 
如 表 2 所 示 。 仿 真 计算 条 件 与 1.3 TIS, US 
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同样 选择 RAE2822 jf 7 38 70 


Casel 
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Fig. 4 Bump Configuration optimization interval selection 
表 2 鼓 包 优化 区 间 参 数 


Table 2 Bump Configuration optimization interval 


pa 


Casel 
Case2 
Case3 


rameters 


X/C 
0.35-0.57 
0.42-0.64 

0.5-0.71 


优化 后 的 升 阻 系数 结果 如 表 3 所 示 , 不 同 Case 


下 的 优化 结果 均 呈 现 不 同 程度 的 减 阻 效果 。 当 翼 型 
e 2 区 间 时 ， 具 有 最 好 的 减 阻 
效果 。 在 升力 系数 变化 不 大 的 情况 下 (增加 了 约 
1.59%)， 阻 力 系 数 下 降 了 约 6.15%， 升 阻 比 提高 


上 表面 变形 范围 取 Cas 


8.25%。 远 高 于 Case 1 


1.3906e-02 


Table 3 
optima 
Cd 
Base 
1.3663e-02 
Casel 
1 1.7596 
1.3051e-02 
Case2 
1 6.15% 
1.3467e-02 
Case3 


13.1695 


| conditions 
CI 
7.0339e-01 
7.0456e-01 
1 0.17% 
7.1457e-01 
1 1.59% 
7.1012e-01 
1 0.96% 


区 间 和 Case 3 [X JR] . 
ak 3 不 同 优化 条 件 下 的 气动 特性 


Aerodynamic characteristics under different 


LD 
50.5805 
51.5691 
1 1.9596 
54.7531 
1 8.25% 
52.73157 
1 4.2596 
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图 5 不 同 优化 条 件 下 的 气动 特性 


Fig.5 Aerodynamic characteristics under different optimal 


conditions 


均 表 现 为 : 由 


根据 压力 系数 分 布 曲 线 ， 图 


6， 所 有 优化 结果 


最 初 的 一 道 强 激 波 变 为 优化 后 的 一 道 


弱 激 波 ， 辟 型 表面 变形 实现 了 对 激 波 的 弱化 。 对 于 


AARRE 


Pg. UAH, Case 1 的 激 波 出 
MMAR ECC ANCK IAEA. 
数 分 布 接近 。 对 于 Case 2 和 Case3， 从 图 中 可 以 看 


压力 系 


出 激 波 明显 减弱 ， 同 时 激 波 出 现 的 位 置 相 对 原始 避 


型 后 移 。 


变形 区 域 中 心 附近 时 ， 


根据 计算 结果 ， 可 以 认为 当 激 波 位 置 处 于 鼓 包 


(有 最 佳 的 减 阻 效果 ， 如 


Case 2 的 计算 结果 ;， 当 激 波 位 置 处 于 鼓 包 变 形 区 域 


中 心 偏向 前 缘 端 时 ， 仍 有 


激 波 位 置 处 于 鼓 包 变形 
阻 效 果 最 小 。 


区 


定 效果 ， 如 Case 3; 当 


ur 


中 心 偏向 后 缘 端 时 ， 减 


(a) Pressure Coefficient 


(b) Detail View 
6 压力 系数 分 布 
Fig.6 Pressure coefficient distribution 
3 不同 气动 状态 下 的 优化 结 
3.1 攻 角 变化 对 气动 特性 的 影响 


自 适应 鼓 包 的 设计 目的 在 于 飞行 器 可 以 根据 自 
身 飞 行 状态 的 不 同 ， 实 时 的 改变 鼓 包 的 构 型 ， 以 保 
证 在 所 需 的 飞行 条 件 下 ， 机 经 始终 具有 较 好 的 气动 
特性 。 为 此 ， 需 要 对 不 同 气 动 状态 下 的 鼓 包 构 型 进 
行 优化 ， 得 到 一 系列 的 鼓 包 形状 曲线 。 

以 变 攻 角 模 式 为 例 ， 为 了 得 到 不 同 气动 状态 下 
的 鼓 包 曲 线 ， 分 别 选取 不 同 的 攻 角 ， 在 优化 区 间 
Case2 区 域 下 进行 优化 分 析 。 优 化 攻 角 状态 分 别 选 
取 AOA=1.8”，AOA=2.31”，AOA=2.8”。 仿真 
计算 条 件 与 2.3 小 结 相 同 ， 经 型 同样 选择 RAE2822 
超 临界 辟 型 。 如 图 7 所 示 ， 鼓 包 高 度 随 优化 攻 角 的 
增 大 而 增 大 ， 最 大 高 度 变 化 小 于 3%。 
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图 7 不 同 攻 角 下 的 鼓 包 优化 构 型 对 比 


Fig.7 Comparison of optimized configuration of bump at 


different angles of attack 
ENAA FAROA ETRE EAE 
表明 , 在 不 同 攻 角 下 鼓 包 均 能 改善 豆 型 的 气动 性 能 ， 
表现 在 阻力 系数 减 小 ， 升 阻 比 增 大 。 同 时 ， 相 对 于 
阻力 系数 、 升 阻 比 的 变化 ， 升 力 系 数 并 没有 出 现 大 


幅 变 化 ， 基 本 保持 稳定 (改变 量 小 于 2% )。 

在 翼 型 攻 角 较 小 时 (AOA=1.8” 0, XU EX 
无 明显 的 激 波 ， 因 此 优化 后 的 鼓 包 较 小 ， 对 原始 避 
型 的 气动 特性 影响 不 大 ， 升 力 系 数 、 阻 力 系数 的 变 
化 小 于 1%。 随 着 攻 角 的 增加 (AOA=2.31”， 
AOA=2.8”), 器 型 上 表面 出 现 较 强 的 激 波 , $x017T 
始 起 到 明显 的 减 阻 效 果 ， 阻 力 系数 分 别 减少 6.15% 
和 10.28%。 同 时 ， 随 着 攻 角 的 增加 ， 鼓 包 的 高 度 也 
随 之 增加 。 

表 4 不 同 优化 条 件 下 的 气动 特性 


Table 4 Aerodynamic characteristics under different 


j 


optimal conditions 


AOA CL CD LD 

2822-1.8 6.10e-01 1.165e-02 52.3664 

6.13e-01 1.158e-02 52.9243 
Bump-1.8 

1 0.5196 } 0.55% 11.0795 
2822-2.31 7.03e-01 1.391e-02 50.5805 

7.15e-01 1.305e-02 54.7531 
Bump-2.31 

t 1.59% } 6.15% 1 8.25% 
2822-2.8 7.90e-01 1.757e-02 44.944 

8.04e-01 1.577e-02 50.984 
Bump-2.8 

11.7896 1 10.2896 1 13.4496 
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如 图 8 BUTEER HR ZERTUUE RH, EIKER, Sx 
f f RAJ. EKER, KAE D en 
RARI. WE 9 所 示 ， 在 大 于 优化 攻 角 的 
条 件 下 ， 对 比 原始 翼 型 ， 增 加 鼓 包 后 的 升 阻 比 有 较 
大 幅度 的 提升 。 同 时 ， 从 图 中 可 以 看 出 ， 对 茶 一 特 
定 攻 角 下 的 优化 结果 ， 当 攻 角 小 于 该 优化 攻 角 时 ， 
气动 特性 表现 较 原 始 翼 型 差 。 
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图 8 不 同 优化 攻 角 下 的 升 阻 极 曲线 
Fig.8 Drag polars for the clean airfoil and different bump 


shapes at the design angles of attack of the bumps 
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图 9 不 同 优化 攻 角 下 升 阻 比 随 攻 角 的 变化 


Fig.9 L/D vs AOA for the clean airfoil and different bump 

shapes at the design angles of attack of the bumps 

考察 多 个 攻 角 下 的 鼓 包 优 化 结果 ， 如 图 10a) 
所 示 ， 利 用 1.2 节 的 优化 方法 ， 对 多 个 攻 角 下 的 自 
EMRO ETRE, FARKE REBAH 
性 。 可 以 看 出 ， 通 过 对 不 同 攻 角 下 的 鼓 包 构 型 进行 
优化 ， 在 大 攻 角 下 的 升力 有 较 大 提升 ， 最 大 升 阻 比 
提高 了 4.4%， 最 大 升 阻 比 出 现 的 位 置 从 1.8” 攻 和 角 
提高 到 了 2.4” 攻 角 。 同 时 相 较 于 干净 的 RAE2822 
如 型 ， 通 过 对 不 同 攻 角 下 的 鼓 包 构 型 进行 优化 ， 使 
得 变形 后 的 到 型 可 以 在 一 个 大 的 攻 角 范围 内 有 更 高 
的 升 阻 比 。 如 图 所 示 ， 优 化 后 的 辟 型 在 1.8”-2.6° 
攻 角 范围 内 ， 升 阻 比 始终 大 于 原始 辟 型 。 图 b 为 对 
应 攻 角 下 优化 后 得 到 的 鼓 包 最 大 变形 高 度 ， 对 比 图 
a 和 图 b, 可 以 看 出 ， 当 攻 角 小 于 1.8” 时 ， 即 小 于 初 
始 豆 型 最 大 升 阻 比 状态 时 ， 鼓 包 变 形 可 以 忽略 ， 对 
气动 特性 基本 没有 影响 。 当 攻 角 大 于 1.8” 时 , 优化 
后 的 鼓 包 最 大 变形 高 度 与 攻 角 近似 呈 线 性 增长 ， 气 
动 特性 也 有 较 大 的 提升 。 
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(a) Lift to drag ratio 
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(b) Optimized results of bump height at different AOA 
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10 不 同 攻 角 下 优化 前 后 的 优化 结果 对 比 


Fig. 10 Result curve before and after the optimization. 


可 见 ， 使 用 能 够 随 飞行 状态 变化 而 变化 的 自 适 
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的 工 况 条 件 下 都 能 具备 最 优 的 气动 特性 。 


3.2 优化 结果 在 非 设 计 点 的 气动 特性 

飞机 在 飞行 过 程 中 ， 由 于 阵风 等 因素 ， 飞 行 状 
态 会 在 短 时 间 内 迅速 变化 ， 需 要 研究 同一 鼓 包 构 型 
在 不 同 飞 行 条 件 下 的 气动 特性 。 图 11 分 别 为 在 不 同 
攻 角 下 优化 得 到 的 辟 型 在 攻 角 1.8" .231^ 28^, 
马赫 数 0.729 下 的 升力 系数 和 升 阻 比 ， 不 同 攻 角 下 


的 优化 结果 利用 鼓 包 高 度 Y/C 代表 。 


如 图 所 示 ， 在 设计 点 状态 的 优化 结果 由 黑 圈 标 
出 ， 其 升力 系数 相对 于 干净 翼 型 变化 不 大 ， 但 具有 
最 高 的 升 阻 比 ， 说 明 阻力 系数 大 幅 减 小 。 当 鼓 包 高 
度 小 于 设计 点 最 优 高 度 时 ， 升 力 系数 基本 不 变 ， 升 
阻 比 相对 于 原始 器 型 有 一 定 的 提高 ， 提 高 幅度 小 于 
设计 点 。 当 鼓 包 高 度 大 于 设计 点 最 优 高 度 时 ， 升 力 
系数 、 升 阻 比 均 迅速 减 小 。 因 此 ， 只 有 当 鼓 包 高 度 


小 于 等 于 最 优 高 度 时 ， 鼓 包 能 够 起 到 减 阻 的 效果 。 
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(a) Lift coefficient VS Height of bump 
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(b) Lift to drag ratio VS Height of bump 
11 气动 性 能 随 鼓 包 高 度 的 变化 规 得 
Fig. 11 Aerodynamic performances of different bump shapes 
考察 攻 角 1.8”、2.31”、2.8” 下 优化 得 到 的 新 
辟 型 ， 在 气动 状态 为 攻 角 2.31”、 马 严 数 0.729 时 
的 气动 特性 , 图 12 a. b 分 别 为 其 压力 系数 分 布 和 
气动 特性 变化 率 。 可 以 看 出 ， 在 攻 角 2.31” 状态 下 
得 到 的 优化 结果 (具有 最 优 的 鼓 包 高度 ) 具有 最 优 
的 气动 特性 ,阻力 系数 降低 了 6.15%; 在 1.8” 攻 角 
下 优化 得 到 的 器 型 〈 鼓 包 高 度 小 于 最 优 高 度 )， 在 
2.31” 攻 角 状 态 下 ， 鼓 包 对 流 场 的 影响 很 小 ， 辟 型 
气动 性 能 的 改善 较 小 ， 阻 力 系数 降低 了 2.51%; 在 
2.8” 攻 角 下 优化 得 到 的 性 型 ( 喜 包 高 度 大 于 最 优 高 
ED, 在 2.31” 攻 角 状 态 下 , 气流 在 鼓 包 前 方 已 经 由 
于 受到 压缩 出 现 第 一 道 激 波 ， 在 鼓 包 之 后 又 出 现 第 
TKA, KEZE TERARI, MIRA 
增加 了 23.196. 


E 


----- SCB OptAOA-1.8? 


S | SCB OptAOA-2.31? 
e SCB OptAOA-2.8? 
-1.5 i 
0 02 04 06 08 1 
X/C 


(a) Pressure coefficient distribution 
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(b) Aerodynamic performances 
图 12 不 同 鼓 包 构 型 在 马赫 数 0.729， 攻 角 2.31" 条 件 下 的 
气动 特性 
Fig. 12 Aerodynamic performances of differentbump shapes 
at a design point of Ma = 0.729 andAOA =2.31° 

因此 ， 为 充分 发 挥 鼓 包 的 减 阻 效 果 ， 同 时 不 影 
响 喜 型 的 升力 特性 , 采用 可 主动 变形 的 自 适应 鼓 包 ， 
根据 需要 实时 改变 自身 构 型 ， 是 一 个 有 效 的 选择 。 


4 激 波 控制 鼓 包 减 阻 机 理 


现代 超 临界 机 豆 在 设计 过 程 中 ， 要 求 辟 型 设计 
点 压力 分 布 为 无 激 波 压力 分 布 。 但 是 ， 由 于 传统 机 
辟 辟 型 不 可 变 ， 无 激 波 特征 只 能 维持 在 一 个 较 罕 的 
飞行 工 况 内 。 随 着 速度 或 者 攻 角 的 增加 ， 机 辟 表 面 
会 出 现 较 强 的 激流， 合适 的 鼓 包 形状 可 以 改变 辟 型 
局 部 构 型 ， 减 小 波 阻 。 根 据 激 波 控制 鼓 包 〈SCB ) 
相关 研究 ， 鼓 包 减 阻 的 原理 主要 有 两 种 形式 : 1) 入 
形 激 波 减 阻 原理 ; 2) 等 烂 压缩 减 阻 原理 。 两 种 原理 
分 别 对 应 不 同 的 鼓 包 构 型 。 前 述 鼓 包 建 模 过 程 中 ， 
通过 对 NURBS 曲线 控制 点 进行 限制 ， 保 证 了 鼓 包 
跟 原 始 辟 型 之 间 光 滑 过 度 ， 使 得 救 包 构 型 为 凹 - 凸 - 
止 构 型 ， 因 此 可 以 利用 等 商 压 缩 原理 进行 解释 。 当 
气流 流 过 连续 光滑 的 鼓 包 曲面 时 ， 气 流连 续 转 折 ， 
产生 无 数 道 微弱 的 压缩 波 ， 这 些 压缩 波 聚 集 在 一 起 
形成 了 一 道 较 弱 的 斜 激 波 呈 。 这 一 过 程 增加 了 辟 型 
上 表面 的 等 商 压缩 范围 ， 使 得 表面 气流 压缩 更 为 充 
分 ， 弱 化 了 激 波 强度 ， 减 小 了 波 阻 。 

如 图 13 所 示 , 随 着 优化 条 件 的 改变 ( 变 攻 角 )， 
鼓 包 减 阻 的 效果 越 来 越 强 ， 主 要 减 小 的 阻力 分 量 为 
激 波 阻力 。 
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图 13 不 同 优化 鼓 包 的 减 阻 效果 
Fig. 13 Components of drag gain for different bump shapes at 


a design point 
UA AOA=2.8” 为 例 ， 从 图 14 速度 云图 可 
以 看 出 ， 使 用 自 适 应 鼓 包 ， 可 以 明显 改善 非 设 计 点 
的 气动 特性 。 优 化 后 的 鼓 包 增 加 了 辟 型 上 表面 的 等 
入 压缩 区 ， 对 上 表面 气流 压缩 的 更 充分 ， 因 此 使 得 
激 波 出 现 位 置 后 移 ， 同 时 强度 减弱 。 激 波 位 置 从 
56%c 位 置 后 移 到 62%c 的 位 置 。 
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(a)Isobars in the shock region ofthe clean airfoil 
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(b) Isobars in the shock region of the airfoil 
with an optised bump 
图 14 激 波 区 域 速 度 分 布 
Fig. 14 Isobars in the shock region 
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5 基于 SMA 的 自 适应 鼓 包 概念 设计 


5.1 自 适应 鼓 包 设计 原理 

形状 记忆 合金 SMA) 具有 独特 的 形状 记忆 效 
应 效应 号 中， 广泛 应 用 于 多 种 工业 领域 。 通 过 专门 
的 时 效 处 理 ， 可 以 得 到 具有 双 程 记忆 效应 的 合金 材 
料 。 对 记忆 合金 结构 施加 不 同 的 控制 温度 ， 结 构 能 
够 在 高 温 构 型 和 低温 构 型 之 间 连 续 变 形 (> 791, 

利用 形状 记忆 合金 的 双 程 形状 记忆 效应 , 可 以 实 
M HEME (Adaptive Bump) 结构 设计 。 如 图 15 
所 示 ， 鼓 包 主动 变形 的 实现 主要 依靠 三 个 方面 : D 
SMA 蒙 皮 构 成 的 鼓 包 ; 22 SMA 加 热 装 置 (电磁 涡 
流 加 热 等 ); 3) SMA 冷却 (空气 冷却 )。 


空气 冷却 


SMAS, 


“SMA Inde 
15 SMA 自 适应 鼓 包工 作 原 理 
Fig.15 Working principle of SMA adaptive bump 
5.2 SM 鼓 包 变形 仿真 
为 研究 SMA 自 适 应 鼓 包 结构 设计 的 可 行 性 ， 对 
SMA 条 带 进行 仿真 研究 ，SMA 本 构 模 型 采用 
Lagoudas 三 维 唯 象 理 论 横 型 请 21。 通 过 SMA 材料 
的 弹性 应 变 、 热 膨胀 应 变 以 及 由 材料 相 变 引起 的 应 
变 来 描述 结构 变形 过 程 中 产生 的 应 变 。 念 真 流程 如 
图 16 所 示 
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图 16 形状 记忆 合金 变形 仿真 程序 流程 图 


Fig.16 Flow chart of simulation program for deformation of 


shape memory alloy 


界 条 

的 预 

YT TR 证 
> 氏 体 


SMA 材料 的 主要 参数 设置 如 下 表 所 示 
表 1 SMA 材料 参数 


Table 1 Shape memory alloy material parameters 
材料 参数 数值 

1 未 变形 条 带 尺 十 140*20*0.5mm 

2 奥 氏 体 杨 氏 模 量 70 GPa 

3 ” 马 氏 体 杨 氏 模 量 30 GPa 

4 ” 正 相 变 开始 温度 Ms 291K 

5 正 相 变 结束 温度 Mf 271K 

6 逆 相 变 开 始 温度 As 295K 

7 逆 相 变 结束 温度 Af 315K 
如 图 17 所 示 ，SMA 条 带 仿真 过 程 中 ， 念 真 边 
件 为 两 端 固 支 ， 在 条 带 长 度 方向 上 有 一 个 -0.01 


应 变 ， 温 度 载荷 加 载 范围 从 290K 到 350K， 以 
条 带 能 够 经 历 一 个 完整 的 相 变 循环 〈 材 料 由 马 
相 到 奥 氏 体 相 变化 )。 


U, Magnitude 
*5.777e-04 


44.044e-04 


*3.466e-04 
-2.888e-04 
*2.311e-04 
*1.733e-04 
*1.155e-04 
T5.777e-05 
+0.000e+00 


(a) 295K 


U, Magnitude 
T*6.299e-03 


*5.669e-03 


43.149e-03 


46.299e-04 
+0.000e +00 


化 的 
长 度 
spit 
变化 


(b) 350K 
图 17 形状 记忆 合金 条 带 有 限 元 仿真 
Fig.17 Finite element simulation of SMA ribbon 
图 18 ARHED REFER Ic E B da 3E 26 
曲线 。 结果 表明 ,SMA 条 带 最 大 挠 度 达 到 条 带 


的 4.4%。 同 时 ， 从 逆 相 变 开始 温度 As (295K) 
相 变 结束 温度 Af (315K) 之 间 , 变形 近似 呈 线 性 
， 具 有 较 好 的 可 控 性 。 
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图 18 形状 记忆 合金 条 带 温度 -位 移 曲 线 
Fig. 18 


SMA ribbon temperature displacement curve 


6 结论 
(1) 使 用 NURBS 曲线 对 自 适 应 鼓 包 进 行 参 


数 化 建 模 。 再 此 基础 上 ， 应 用 isight 优化 平台 集成 
了 Matlab 鼓 包 曲 线 生 成 、Gamibt 网 格 更 新 、Fluent 


流动 求解 模块 ， 搭 建 了 用 于 自 适 应 鼓 包 优 化 的 求解 
2 


可 以 利用 遗传 算法 进行 优化 求解 。 
C20 利用 鼓 包 优化 平台 研究 了 RAE2822 WA 


不 同 喜 包 变 形 位置 对 减 阻 的 影响 ， 确 定 了 最 优 的 变 


形 区 域 。 研 究 表 明 ， 当 激 波 位 置 位 于 鼓 包 中 心 时 具 
有 更 好 的 减 阻 效果 。 对 于 RAE2822 A, 42%-65% 
弦 长 处 为 自 适应 鼓 包 最 优 的 位 置 。 


最 优 豆包 构 型 。 结 


(3) 优化 得 到 了 RAE2822 AWER EAA F 
表明 ， 随 着 攻 角 的 增加 ， 鼓 包 


减 阻 的 效果 越 来 越 明 显 。 在 攻 和 角 AOA=2.31” 时 ， 
IRBE 6.15%; 攻 角 AOA=2.8° 时 ， 减 阻 10.28%。 


AR 


3H 
包 


适应 鼓 包 设计 方案 。 利 | 
对 条 带 鼓 包 变形 ; 


(4) 为 充分 发 挥 鼓 包 的 减 阻 效果 , 同时 不 影响 


型 的 升力 特性 ， 需 要 采用 可 主动 变形 的 自 适 应 鼓 
,根据 气动 要 求实 时 改变 自身 构 型 。 


(5) 提出 了 一 种 基于 形状 记忆 合金 材料 的 自 
三 维 SMA 本 构 材 料 程 序 


行 了 仿真 ， 鼓 包 最 大 变形 量 达到 


性 关系 ， 理 论 上 可 以 满足 自 适应 鼓 包 的 变形 要 求 。 
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